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ABSTRACT. Magnetic properties of a ternary-spin Ising model on the decorated square lattice are studied 
within a generalized decoration-iteration transformation. Depending on the mutual ratio between exchange inter- 
actions and the single-ion anisotropy, there apprear six different phases in the ground state. The magnetic order 
of these phases together with the critical behaviour and corresponding magnetization curves are discussed in detail. 



INTRODUCTION 

Decorated ferrimagnetic Ising models have enjoyed 
an immense research interest for many years. Besides 
exactly solvable two-sublattice Ising model [T], much 
effort has been recently devoted to the decorated Ising 
square lattice composed of three magnetic atoms of 
different magnitudes [2]. However, it is noteworthy that 
only the compensation phenomenon of this model has 
been exactly examined yet, while the critical behaviour 
and thermodynamic quantities have been calculated 
with Bethe-Peierls approximation [2]. Owing to this 
fact, the aim of this work is to show a simple way for 
obtaining complete exact results for this system. 

MODEL AND METHOD 

In this article, we shall consider a ternary-spin Ising 
model on the decorated square lattice (see Fig. [I), 
where the lattice sites labeled by white circles are occu- 
pied by A atoms with the spin Sa = 1/2 and the dec- 
orating sites denoted by black and hatched circles are 
occupied by B and C atoms with the spins Sb = 3/2 
and Sc — 5/2, respectively. The total Hamiltonial of 
the system reads: 

H = JrmS:S* m +J2 J m S Z t S Z n -DY,{Stf , (1) 

where the first two summations are carried out over the 
nearest-neighbour A— B and A— C pairs, respectively, 
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S B = 3/2 



S c = 5/2 




and the last one runs over the decorating B and C 
sites. In Eq. S? and S% (k — m,n) denote the 
Ising spin variables of A atoms and B or C atoms de- 
pending on whether k is selected as k = m or k = n. 
Next, Jik stands for the exchange interaction between 
the nearest A— B neighbours if k = m and the nearest 
A— C neighbours if k = n. Finally, the parameter D 
stands for the single-ion anisotropy acting on B and C 
atoms. To proceed further, we introduce the general- 
ized decoration-iteration transformation [3] : 

s k 

exp (f3Dn 2 ) cosh [(3nJ lk (£? + S* +1 )] 

n— — Sk 

= P a exp(j3R a SfS? +1 ) , for a = h,v (2) 

where /3 = l/(k-QT) and the subscript a specifies the 
bond of the lattice to which the transformation @ is 
applied. Namely, Eq. ([2]) is used for horizontal bonds 
if k = m and then one puts a = h, while for k = n it is 
applied for vertical bonds and thus a = v. To simplify 
further calculation, we shall assume that interactions 
Ji m and Ji n take exclusively positive values Jab and 
J ac, respectively. By the use of Eq. ([2]), it is then quite 
straightforward to derive the equality: 

Z(J3, Jab, Jac,D) = P N Z {f3, R hl R v ) , (3) 

which relates the partition function Z of the consid- 
ered model and that one Zq of the original anisotropic 
spin- 1/2 Ising square lattice with nearest-neighbour 
couplings Rh and R v in horizontal and vertical direc- 
tion, respectively. In above, N denotes the total num- 
ber of A atoms. From the Eq. ((3| one can readily 
calculate all important quantities utilizing the stan- 
dard thermodynamic relations. In the present arti- 
cle, we shall turn our attention to the analysis of the 
ground state and the critical behaviour. In addition, it 
is also useful to investigate the sublattice magnetiza- 
tion m A = (Sf ) , m B = (S m ) , m c = {S n }, where (. . .) 
stands for the standard canonical average, and the to- 
tal magnetization reduced per one atom of the original 
Ising model m = M/N — vtla + tob + mc, as well. Fi- 
nally, let us mention a critical condition, which follows 
from Onsager's exact solution [4]: 



Fig. 1. Part of the decorated Ising square lattice. The 
white circles denote A atoms, while the black and hatched 
ones represent decorating B and C atoms, respectively. 



smh{p c R h /2) smh{p c R v /2) = 1, 



(4) 



where (3 C = l/(fc B T c ) and T c denotes the critical tem- 
perature of the model system. 



RESULTS AND DISCUSSION 



Now, we proceed to a brief discussion of the most 
interesting numerical results obtained for the ground- 
state and finite-temperature phase diagram. Directly 
from Fig. [2J where the spin order drawn in square 
brackets shows a typical spin configuration [5a, 5b, 5c] 
in the ground state, one easily finds that A atoms re- 
main in the spin state 5a = 1/2 in the whole pa- 
rameter space, while the spin order of decorating B 
and C atoms is basically modified by the single-ion 
anisotropy D and the ratio Jac/Jab- Actually, a first- 
order transition D* = — 0.5Jab associated with the 
spin crossover of B atoms 5b = —1/2 <-» —3/2 oc- 
curs regardless of Jac/Jab as D reinforces. On the 
other hand, the Jac / Jab strengthening clearly causes 
spin changes of C atoms from the lowest 5c = —1/2 
towards the highest 5c = —5/2 available spin state. 
Therefore, two further first-order transitions J AC = 
—2D and J A ^ = — 4D corresponding to spin crossovers 
5 C = -1/2 <-> -3/2 and 5 C = -3/2 -5/2, respec- 
tively, can be found within the D < 0.0 region. In 
consequence of this, there appear in total six different 
ferrimagnetically ordered phases in the ground state 
depending on the ratio between D, Jab and Jac- 

Next, let us make some comments on the finite- 
temperature phase diagram displayed in Fig. O which 
shows the variation of critical temperature with the 
single-ion anisotropy D for several values of the ra- 
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Fig. 2. The ground-state phase diagram. 




Fig. 4. The thermal variations of the total magnetization 
for Jac/Jab =1.5 and several values of D. 

tio Jac/Jab- As one can see, the critical tempera- 
ture monotonically decreases from its maximum value 
at D — > 00, and tends towards its global minimum 
achieved in the limit D — > — 00 regardless of the ra- 
tio Jac/Jab- Thus, the spin system remains long- 
range ordered in the ground state independently of 
the Jac / Jab strength in agreement with the previous 
ground-state analysis. Open, crossed and black circles 
denote in Fig. [3] those critical temperatures at which 
spin crossovers 5c = —1/2 <-» —3/2, 5c = —3/2 <-» 
—5/2 and 5b = — 1/2 <-> —3/2, respectively, emerge in 
the ground state. Besides, symbols arising from mix- 
ing of relevant circles represent special critical points 
that correspond to the coexistence of four different fer- 
rimagnetic phases (compare Figs. [2] and [3]). 

In Fig. [4[ we show some typical thermal variations 
of the total magnetization for Jac/Jab = 1-5 and 
several values of D. As one can see, the plotted re- 
sults are completely consistent with those presented in 
Figs. [2] and [H Indeed, the initial value of |m| takes 
one of four possible values |m| = 1/2,3/2,5/2 or 7/2 
corresponding to [±, -|, -±], [|, -§], [|, -§, -§] 
or [i — |, — |] phases, respectively. Furthermore, there 
are also three special cases starting at |m| = 1,2 and 3 
(labeled by dashed lines in Fig. [4j) , which correspond 
to two-phase coexistence at T = 0. 

CONCLUSIONS 
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Fig. 3. The critical temperature versus the sigle-ion 
anisotropy D for several values of the ratio Jac / Jab ■ 



In conclusion, it should be mentioned that the 
model under investigation shows rich critical behaviour 
and many other interesting features, which require 
further detailed analysis. Another interesting results 
concerning this model will be therefore presented in 
the more comprehensive review to be published in the 
near future. 
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